Anomalous modes drive vortex dynamics in confined Bose-Einstein condensates 
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The dynamics of vortices in trapped Bose-Einstein condensates are investigated both analytically 
and numerically In axially symmetric traps, the critical rotation frequency for the metastability 
of an isolated vortex coincides with the largest vortex precession frequency (or anomalous mode) 
in the Bogoliubov excitation spectrum. As the condensate becomes more elongated, the number of 
anomalous modes increases. The largest frequency of these modes exceeds both the thermodynamic 
critical frequency and the nucleation frequency at which vortices are created dynamically. Thus, 
anomalous modes describe not only the critical rotation frequency for creation of the first vortex in 
an elongated condensate but also the vortex precession in a single-component spherical condensate. 
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Bose-Einstein condensation (BEC) and superfluidity 
are two entangled core issues of low-temperature physics, 
and recent experimental developments now allow us to 
study them in a nearly ideal system: a dilute gas of al- 
kali atoms with well-understood interactions [|l|-^ . Many 
key aspects of BEC have been clarified since the first ex- 
perimental observations, and during the past year much 
attention has been given to manifestations of superflu- 
idity. The 'scissors mode' of excitation of trapped con- 
densates Q implies the irrotational flow characteristic of 
superfluids. Quantized vortices, which have long been 
known as fundamental excitations of superfluid helium 
and superconductors, have been observed directly in one- 
lH,|j and two-component trapped Bose gases. 

Multiple vortices have recently been generated in con- 
flned single-component condensates B by rotating a 
weakly anisotropic trap at an angular frequency Q. This 
approach resembles the classic rotating-bucket experi- 
ments on liquid helium Ig] , where vortices first appear |9|] 
at a critical frequency f2c above which the presence of 
the vortex lowers the total free energy [[lOyil| ] . There is a 
discrepancy, however, between the observations of Madi- 
son et al. [g and simple theoretical considerations based 
on the equilibrium solution of the Gross-Pitaevskii (GP) 
equation, which has otherwise been remarkably success- 
ful in describing condensate behavior ||l|; in particular, 
the frequency at which vortices are first observed is sig- 
nificantly higher than equilibrium estimates of Qc- 

The resolution of this discrepancy is the subject of the 
present paper. In particular, we find that the 'anoma- 
lous modes' of the Bogoliubov spectrum determine the 
critical rotation frequencies at which vortex arrays are 
observed in the ENS experiment M ; they also determine 
the frequency of precession of the vortex core observed 
at JILA §. 

In trapped condensates, several factors influence the 
critical rotation frequency for appearance and stability 
of vortices, and for their subsequent dynamical motion: 



(1) An energy barrier at the surface of the conden- 
sate [|l2| inhibits vortex penetration of its interior. Thus, 
vortices are predicted to nucleate spontaneously at a 
frequency n^, > f2c ||l3| , whose value coincides with 
the criterion for the Landau instability of surface exci- 
tations ]14|-|l^ and with the effective disappearance of 
the surface barrier [|l^. We show below that, while Qu 
exceeds ilc for the cigar geometry of the ENS experi- 
ments B, it is still lower than the frequency at which 
vortices are first observed. 

(2) Stringari pj] has shown that Qc increases with T, 
reaching a maximum close to the BEC transition T^. The 
temperature T ~ O.STc required to match the observed 
critical frequency, however, is much larger than the tem- 
perature at which the experiments are performed. 

(3) A condensate with a vortex first becomes sta- 
ble against small perturbations at the metastability fre- 
quency flm- One definition for metastability is that the 
energy per particle must not decrease under infinitesimal 
displacements of the vortex from the condensate center; 
for a disk-shaped condensate, this yields O^ = fi^c in 
the Thomas- Fermi (TF) limit ^. Equivalently, the Bo- 
goliubov excitation spectrum must contain only modes 
with positive energy. A non-rotating axisymmetric con- 
densate with a singly quantized vortex has at least one 
"anomalous mode" ; such modes have a negative excita- 
tion energy Sa and a positive norm |lj,|l^ . For a rotating 
axisymmetric trap, these anomalous modes are Doppler- 
shifted upward by f2, with the metastability frequency 
rirn = niax(|ea|)/?i. As shown below, r2,„ can exceed 
both flc and Q^, in cigar traps. 

This criterion for the onset of linear stability agrees 
well with the large frequency at which vortices first ap- 
pear experimentally in cigar-shaped condensates [pi. It 
also explains why the observed critical frequency is inde- 
pendent of whether the condensate is first cooled and sub- 
sequently rotated, or vice versa. Otherwise, the critical 
rotation should be greater for the cool-and-then-rotate 



scenario, since vortices must surmount the energy barrier 
at the condensate surface, penetrating the cloud only at 
the nucleation frequency Q„ > He- 

(4) The anomalous modes correspond to the precession 
of the vortex core about the origin of the trap |13,M . For 
il < flrn , the vortex moves in the same direction as the 
superfluid flow around the core. Hence, the anomalous 
mode(s) should describe the motion of one-component 
vortices in recent JILA experiments |^,0]- 

The present work links anomalous modes with these 
various observations. We proceed from the zero- 
temperature description of the order parameter (or wave 
function) for a condensate of ^^Rb atoms, as a solution 
time-dependent GP equation 



idt^Pir, t) = {T + V^trap + Vk - ni,) ^(r, t). 



(1) 



with the kinetic energy operator given by T = — ^V^, 
the Hartree field by T4j = 47r77|-0p, and trap potential 
by Vtrap = 5 [A2(1 -I- e,,)x^ + X^il + ey)y^ + z^] , where 
A = u)p/(jjz is the trap anisotropy and e^, ey describe 
(small) departures of the trap from axial symmetry. The 
centrifugal term — ilL^ = ifl (xdy — ydx) appears for sys- 
tems rotating about the z axis at a frequency fi. Energy, 
length, and time are respectively given in harmonic os- 
cillator units hiOz, dz = a/S/Mwj,, and wj^, where Uz is 
the axial trap frequency and M is the atomic mass. Nor- 
malizing the condensate wave function J(ir|'i/;(r, t)|^ = 1 
yields the scaling parameter 77 — Na/dz, where a = 
lOOao « 5.29 nm is the ®'^Rb scattering length |^ and N 
is the number of condensate atoms. As parameters typi- 
cal of the recent experiments at ENS Q and JILA [^,0 , 
we take iV = 1.4 x 10^ with {h'p,iyz) = {ujp,LL!z)/2n — 
(169, 11.7)Hzand(e^,ej,) = (0.03,0.09), and A^ = 3x10^ 
with (7.8, 7.8) Hz and ex = ^y = 0, respectively. 

The stationary solutions of the GP equation, defined 
as local minima of the free energy, are found numerically 
by norm-preserving imaginary time propagation using 
an adaptive stepsize Runge-Kutta integrator. The spa- 
tial dependence of the complex condensate wave function 
employs a discrete-variable representation |g^ based on 
Gaussian quadrature, and is assumed to be even under 
inversion of z. The numerical techniques are described 
in greater detail elsewhere ||l^,|2^. The initial conden- 
sate amplitude is taken to be the TF wave function 
■0TF = y^iiI^^7^^\^^^^^)JI^, where /xtf = ^(ISA^r;)^/^ 
is the chemical potential. The GP equation for a given 
value of Q and N is propagated in imaginary time until 
the fluctuations in both the chemical potential and the 
norm become smaller than 10~^^. 

Once the system reaches equilibrium, its response 
to small disturbances is found by substituting ^ -^ 



-ifj,t 



{^|;- 



ue 



■ ve 



*) in Eq. 



The Bogoliubov 



spectrum of eigenvalues e is found by linearizing in the 
quasiparticle amplitudes u(r) and v{r). 

Figure H shows the critical frequencies Qc/^n and 
rij,/27r as a function of trap anisotropy A. We ignore 



the small in-plane trap distortion (setting ex = ey ~ 0) 
in order to simplify the computation. In this axisym- 
metric system, the thermodynamic critical frequency flc 
is the energy difference between condensates with and 
without a vortex at the trap center (divided by h), since 
the rotation Doppler shifts the former by exactly — fi; 
for the parameters of the ENS experiment, A « 14.44, 
we obtain Slc/27r « 0.4i'p « 68 Hz, in good agree- 
ment with the TF estimate but much lower than the 
observed value riobs/27r « OJi/p « 120 Hz. We find 
that inclusion of the small in-plane anisotropy does not 
appreciably increase the value of fie, in contrast with 
results reported recently |Q. The critical nucleation 
frequency fi^ = min^Snm/^fn) defines the rotation fre- 
quency at which the first Bogoliubov excitation of the 
vortex- free condensate becomes unstable iflil . Wc obtain 



rii./27r « 0.52i^p 



Hz, which is larger than Qc/'2t^ but 



still too low to account for the experimental result. 
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FIG. 1. The numerical values for the thermodynamic crit- 
ical frequency Q,c/2ti (solid curve), the nucleation frequency 
r2^/27r (dashed curve), and the even z-parity anomalous mode 
frequencies (circles) are shown as a function of trap anisotropy 
A — Vp/vz. Parameters are N = lA x W\ u^ — 11.7 Hz, 
and ex = ey = 0. The dotted vertical line corresponds to 
anisotropy relevant to the ENS experiments, and the '®' la- 
bels the frequency at which vortices are first observed. 

Also shown in Fig. ^ are the frequencies of the anoma- 
lous modes |a;a|/27r = |ea|/^ with even z parity (plot- 
ted as positive values). These are Bogoliubov excita- 
tions of the vortex state with negative energy but positive 
norm, and relative angular momentum —fi. For spherical 
or pancake geometries, the spectrum contains only one 
anomalous mode. It describes the precession of the vor- 
tex about the center of the trap, in the same sense as the 
circulation about the core [ po| at a frequency Wq = |f2c 
in the TF limit [|l2|. As A increases and the vortex line 
stretches, however, additional anomalous modes appear 



in the energy spectrum |^,^ , corresponding to the pre- 
cession of a spatiahy deformed vortex hne. Physically, 
flc involves only a straight vortex (note, however, that 
Ref. [g3| finds a deformed vortex for the ground state in 
very elongated condensates), but f2m involves linearized 
deformations. For an elongated condensate (A ;» 1), the 
large curvature of the condensate surface readily induces 
distortions. The Bogoliubov amplitudes are highly local- 
ized radially in the vortex core and oscillate as a function 
of z, with maximum amplitude near the condensate sur- 
face. In the TF limit, n anomalous modes appear above 
a critical anisotropy A„ > yjn{n + l)/2 [p5| . 

Excitation of the anomalous modes can lower the sys- 
tem's free energy, destabilizing the vortex state. Metasta- 
bility of the vortex is guaranteed for rotation frequencies 
exceeding the largest anomalous mode, Vtm = max(|tJa[). 
In the TF limit, 0^ exceeds Vt^ when A > 2 p5[ , 
which is close to the value of A « 2.5 found numer- 
ically (see Fig. |l|). For the ENS parameters, numer- 
ical analysis finds the metastability critical frequency 
ri„i/27r ~ Q.livp K, 124 Hz, very close to the observed 
value riobs/27r « 120 Hz. 

For the spherical trap relevant to the JILA experi- 
ments, the anomalous mode frequency in the TF limit 
has the predicted form \iOa\/uj = |(^/i?) ln(1.96i?/0, if 
we neglect terms of relative order (^/i?) ln(i?/^) « 0.1, 
where R = {IhvjY/^ is the TF radius and £, = 1/R is 
the dimcnsionless healing length. We find \uJa\/2Tr = 
1.58 ± 0.16 Hz, in reasonable agreement with the ob- 
served precession frequency of 1.8 ± 0.1 Hz. The value of 
\uja\j however, is sensitive to the number of condensate 
atoms (ex N~'^/^) and the displacement of the vortex 
from the trap center p3| (ex [1 — {p/RY\~^). Interest- 
ingly, the spectrum computed numerically also contains 
a cownter-precessing (i.e. non- anomalous) mode with fre- 
quency 3.63 Hz and two nodes along z; this may corre- 
spond to the distorted 'rogue' vortices observed in two- 
component systems Q. 

To make closer contact between vortex precession and 
anomalous Bogoliubov modes, the dynamics of a trapped 
vortex are explored by real-time propagation of the GP 
equation (^. Once the ground state has been obtained, 
an off-center vortex with counterclockwise circulation is 
imposed on the condensate wave function at t = by the 
method of phase imprinting |2^ . The vortex is displaced 
by 1.57(ip along x from the trap origin, corresponding to 
1.3 (Um and 6.1 /im for the ENS and JILA condensates, 
respectively. In both cases, the vortices undergo coun- 
terclockwise precession, which is the same sense as the 
circulation around the vortex core. 

For the JILA parameters, the period is found to be 
623 ± 3 ms (the uncertainty reflects the coarseness of 
the spatial and temporal grids), yielding a precession 
frequency of uOp/2Ti — 1.61 ± 0.01 Hz in excellent agree- 
ment with the value of the predicted anomalous-mode 
frequency discussed above. The phonons that are also 



produced by phase imprinting ||2^ only weakly affect the 
vortex motion, and they rapidly decay into unobserv- 
able high-frequency collective modes. The vortex itself is 
slightly curved, as discussed in Ref. [^. No noticeable 
variation in the vortex displacement from the trap center 
(spiraling) was found after three full precessions. 

As shown in Fig. 0, precession of an off-center vortex in 
the ENS trap is associated with a pronounced curvature 
of the vortex line. The precession appears to be most 
rapid near the surface of the condensate; by 12 ms, the 
ends of the vortex have returned to their initial location, 
while the center lags behind by almost 180°. Defining 
the precession by the motion near the surface, one ob- 
tains a frequency of approximately 85 Hz. This value 
(and to some extent the shape of the vortex) is consistent 
with the second-largest anomalous mode in the excitation 
spectrum shown in Fig. |l|, |ijJa|/27r = Q.b\Vp « 87 Hz, 
which has two nodes along z. The excitation of this 
precession mode probably arises from imprinting a cir- 
culation pattern aligned with the z axis on a condensate 
that is highly non-uniform radially; the initial straight 
vortex would have some overlap with all of the 11 even 
z-parity negative-energy modes found above for this con- 
figuration. 
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FIG. 2. Numerical results for the precession of a vortex in 
the ENS trap. Frames in raster order correspond to 4 ms 
increments after the initial phase imprint. Squares, which are 
8dp ~ 6.7 /im on a side, show the integrated density down z; 
the corresponding vortex line (viewed 9° off the z axis toward 
—y) is rendered above each frame. 

The critical frequency for vortex nucleation in elon- 



gated traps is confirmed by numerical solution of the 
time-dependent GP equation. The initially axisymmetric 
ENS condensate is simultaneously distorted {£x,£y 7^ 0) 
and rotated at t — 0. The GP equation is propagated in 
the rotating frame for 500 ms with cither il ~ Q.7u!p or 
0.8u)p. No vortices are produced for the ft = 0.7ujp case. 
As shown in Fig. 0, however, when f2 = O.Sujp vortices ap- 
pear at the condensate surface hy t — 200 ms, and fully 
penetrate the cloud hy t = 300 ms. In the absence of 
dissipation, the vortex motion remains extremely turbu- 
lent; the dynamics of nucleating vortices is discussed in 
greater detail elsewhere ||2^. The ENS condensates are 
held in the rotating trap for 500 ms, presumably allowing 
the vortex structures to approach equilibrium. The sim- 
ulations indicate that dissipative processes not included 
in the present study, such as the interaction of the con- 
densate with the thermal cloud |29[| , are important in the 
formation of the vortex arrays observed experimentally. 




FIG. 3. The ENS condensate rotating ai fl — O.Sup as a 
function of time. Images in raster order correspond to 150 ms 
through 500 ms in 50 ms increments, and are each 10 x 10 /im. 
The condensate density (integrated down z) is shown in the 
rotating frame. 



mode uJa{^) becomes imaginary in the range \uJa\ — 6 '^ 
fl < jwal -f (5, where 26 = \ex — e-y\ [^. The onset of 
metastability in such traps occurs only when fl exceeds 
l^al by the appropriate amount, which may help explain 
the relevance of flm ^ \uJa\ for the ENS experiments. 

In summary, anomalous modes are interpreted as defin- 
ing both the vortex precession in the JILA experi- 
ments [g| and the critical frequency for the appearance 
of the first vortex in the ENS experiments H . 
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